Pauli-Limited Superconductivity in Small Grains 
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We report on an exploration of the mean-field phase diagram for Pauli-limited superconductivity in 
small metallic grains. Emphasis is placed on the crossover from the ultra-small grain limit where 
superconductivity disappears to the bulk thin-film limit as the single-particle level spacing in the 
grain decreases. We find that the maximum Zeeman coupling strength compatible with supercon- 
ductivity increases with decreasing grain size, in spite of a monotonically decreasing condensation 
energy per unit volume. 
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I. INTRODUCTION 

Mesoscopic physics may be broadly defined as the 
study of phenomena which depend fundamentally on the 
finite-size of a system, even when that size substantially 
exceeds characteristic length scales associated with mi- 
croscopic degtees of freedom. By this definition, recent 
experimentsoQ in which strong parity effects are seen 
in superconducting islands containing ~ 10 9 electrons 
highlight the robustness of superconductor mesoscopics; 
pairing physicsEm causes observable differences between 
metallic grains with 10 9 electrons and grains with 10 9 + 1 
electrons. Recent progress has enabled studieSj-pLmuch 
smaller systems. Ralph, Black, and TinkhamETO have 
demonstrated that it is possible to make single-electron 
transistors with superconducting islands that are only a 
few nanometers in radius. These systems contain only ~ 
10 4 to 10 5 electrons and have a mean energy level spacing 
6 which can be larger than or comparable to Ao, the zero- 
temperature superconducting gap fa bulk samples. As 
early as in 1959 Anderson observed!?] that superconduc- 
tivity cannot occur in small grains when the limit S ~ Ao 
is reached. Experimental realization of such ultra-small 
systems has opened the physics of superconductivity in 
this regime to experimental-study and inspired substan- 
tial theoretical interest .E3il3 Many aspects of the exper- 
iments can be qualitatively understood using the sim- 
plest possible BCS model of an ultra-small grain in which 
the levels are assumed to be equally spaced and pairing 
is assumed to occur only between identical orbitals.ta 
In this paper we address the influence of Zeeman cou- 
pling on superconductivity in such a model, emphasiz- 
ing the crossover between the ultra-small grain regime 
and the bulk limit where the Chandrasekhar-Ciogston 
paramagnetic limit, Z < Zc = Ao/y/2, applies .U Here 
Z = gfigB/2 is the Zeeman coupling strength, /is is the 
Bohr magneton, and B is the magnetic induction. 

For a constant level spacing spectrum, the single par- 
ticle energy levels measured from the Fermi energy are 
= [n — a)5. Here n = 0, ±1, ±2,. . . , and parity de- 
pendence appears in the quantity a which has the value 



if the number of electrons N is odd and 1/2 if N is 
even, corresponding respectively to chemical potentials 
pinned at and half-way between energy levels. The gap 
equation for a model in which pairing interactions oc- 
cur only between identi«al|-Q|rbitals differs from its BCS 
theory counterpart onlyt3ll3 in the discreteness of the 
quasiparticle level spectrum: 
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f(E n + Z)- f(E n - Z) 



(1) 



where E n = y/ A 2 + i \ j s the dimensionless coupling 
constant and A is determined by solving these equations. 
The upper limit on this discrete sum comes from the en- 
ergy cutoff used in BCS theory to model retarded at- 
tractive interactions and can be expressed in terms of 
Ao using the bulk solution of the zero temperature gap 
equation 



A = 2M<Sexp(-l/A). 



(2) 



Since, at available fields, the magnetic flux through an 
ultra-small grains will typically be much smaller than. 
$0 = hc/2e, coupling to orbital degrees of freedomE£l 
can normally be ignored. The electron spins still couple 
to the magnetic field B, however, splitting the single- 
particle energies, — > (n — a)5 ± Z. 

For small superconducting particles it is essential that 
the occupation probabilities, /, in Eq. (Q) be calculated 
in ensembles including states with only even or odd num- 
bers of particles.B These differ from Fermi occupation 
probabilities only for levels close to the chemical potential 
and only for temperatures ksT ^ 5; at T=0 the even re- 
striction has no effect and the odd restriction has only the 
effect of removing the orbital at the Fermi energy, which 
cannot be paired, from the gap emiation. This model 
was first studied by von Delft et al.E3 to calculate the de- 
pendence of A(T) on 5. They found that A(T) remains 
close to its bulk value until S is close to a critical value 



5 C which is parity dependent: <5° dd /A = 5 
^ ven /A = 2e^ ~ 3.56. Here 7 = 0.577215 



7 ~ 0.89 and 
is Euler's 
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constant. As S — ► S c from below, the critical temperature 
and the zero-temperature gap bath tend to zero. Later 
Braun et al.E-3 and Balian et al.E-3 extended this work to 
the case of finite Zeeman coupling. Ref. |lj concentrated 
on comparison between theory and the experiments by 
Ralph et al.Q, finding good qualitative agreement. Ref. [L5| 
concentrated on the influence of using parity dependent 
distribution functions at finite temperature, predicting 
significant qualitative effects for such quantities as the 
superconducting gap and the magnetization. Here we 
explore the full T — Z phase diagram, examining depen- 
dences on parity and S and focussing on the evolution of 
the T — Z phase diagram toward its bulk limit result as 
8 decreases. 



II. PHASE DIAGRAMS 

The fundamental equation on which our calculations 
are based is the following coupling-constant-integration 
expression for the free-energy .difference between super- 
conducting and normal states:Ej 

- A2 

r^ s n ~\$ 

_ 2 f r A dAIA A-f(E' n+ z)- f (E^z) (3) 

n=1 Jo E n 

In Eq. (g) E' n = + A' 2 . Since, in the second term 
on the right hand side of Eq. (g), A appears only in 
the upper limit of the coupling-constant integration, it is 
easy to verify that the gap equation, Eq. (0), is satis- 
fied at extrema of the condensation energy 
sufficient condition for superconductivity in this model 
is that 9a 2 (^s — ^n)|A 2 =o be negative; this derivative 
changes sign along the surface in (Z, T, 8) space where 
the linearized (A 2 — > 0) gap equation is satisfied. When 
the phase transition is continuous, this surface, ^(T, 8), 
is the boundary of the superconducting region. If, on the 
other hand, (fl a — Q n ) |a 2 =o is positive the grain may 
still be in a superconducting phase if f2 s — VL n is negative 
at a finite value of A. In that case the phase boundary 
is of first order and the corresponding critical value of Z 
is denoted Zi(T,S). 

Below we determine the functions Z 2 (T, 8) and 
Z\ (T, 8) , and thereby show how the phase diagram for 
bulk systems is generalized to systems with finite level 
spacing. 

A. Normal state instability at T = 

It is instructive to start by considering Z 2 (T — 0,(5) 
which can be evaluated analytically by solving the lin- 
earized gap equation |f]). To that end we note that at 
T = the factor [1 - /(£„ + Z) - /(£„ - Z)\ in Eq. ® 
equals zero if £„ — Z is negative and otherwise equals 



one. Hence, as Z is increased an additional pair of states 
(n \,n J,), is blocked from pairing every time a + Z/8 
passes through an integer value. It follows that 

Z 2 (T = 0,S) = (m-a)5 (4) 

where m is the largest integer for which 

M 

22 - — - =ip(M + l-a) - ip(m-a) > -. (5) 

n—m 

Here ip{ x ) is Euler's ip function. Using ip(x) ~ \n(x) for 
large arguments and replacing M using Eq. (|^) we find 
that m is the largest integer for which 

5/A < -exp[-tp(m- a)]. (6) 

The resulting Z 2 (T = 0, 8) is plotted in Fig. [j] and Fig. 
U for evenj-and odd N, respectively. In both cases the 
bulk valued Z 2 {T = 0,8 -> 0) = 1/2 is recovered as 
can be verified by letting m become large in Eqs. (^) 
and (||). Rather than decresing steadily, Z 2 (T — 0,8) 
oscillates around its <5 — > limit with increasing 8. 
The maximum values of Z 2 (T = 0, 8) actually occur for 
8 -> S c (m = 1); Z^ nax /A = e 1 ~ 1.78 for even N 
and Z? ax /A = ei/2 ~ 0.89 for odd N. Both of these 
values exceed Zq/Aq = l/v2) an( i this analytic result 
thus establishes that the Chandrasekhar-Clogston limitEJ 
can be exceeded in the ultra-small particle limit and, as 
we see below, also over a broad range of small particle 
sizes. The decrease of A and of mean-field-theory critical 
temperatures with particle size is not accompanied by a 
corresponding decrease in the maximum allowed Zeeman 
coupling strength. 

B. First-order Transition Phase Boundary at T = 

When A is finite, the pair-breaking condition (Z > 
+ A 2 ) is not satisfied until larger values of Z are 
reached compared to the A = case. Hence, at suffi- 
ciently low temperatures, states with finite A are favored 
over states with A — > 0, causing the superconductor- 
normal transition to be of first order. This physics is 
much the same at finite 5 and in the 5 — ► bulk thin film 
limit. 

We consider first the ultra-small grain limit. At T = 
the integral in Eq. (||) can be evaluated analytically. 

A 2 ™ 

Q s - O n = — - 2 y {max(£' w , Z) - max(£„,Z)} 

n— 1 

~^ln(*/*c) + ^C(3,a)+2(Z-Z a ). (7) 

The first form for the right hand side of Eq. (Q) is exact 
whereas the second form only applies in the ultra-small 
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grain regime: £ 2 > Z > £1, E\ > Z and A/5 <C 1. In 
Eq. (0), C is Riemann's Zeta Function; C(3,0) ~ 1.2021 
and C(3, 1/2) ~ 8.4144. Minimizing H s - Q n we find that 
A 2 = 25 2 ln{5 c /5)/({3,a) and 

Zl (T = 0,5) = Z 2 (T = 0,5) + 6 [1 fffJ (8) 

2C(3,a) 

for <5 — > <5 C . This expression is plotted together with 
exact numerical evaluations of Z\(T = 0,5) in Figs. |l| 
and |2| for even and odd N respectively. Z\{T = 0, 5) and 
Z 2 (T = 0, 5) become equal only as 5 — > 5 C . 

We notice in Fig. [I] that the exact numerical evalua- 
tion recovers the familiar Chandrasekhar-Clogston result 
in the bulk limit: Z X {T = 0,5 -»• 0) = l/\/2A is re- 
covered for both even iV and odd N. The exact numer- 
ical results for Z\(T = 0, 5) were obtained by evaluating 
fl s — fl n as a function of A and locating zero-crossings 
of its minimum as Z is increased. In Fig. ||| we show 
the dependence of (f2 s — Cl n )/N on A for four different 
values of 5 at (T = 0, Z = 0.65Ao) for an even number of 
electrons. For 5/A = 1.4 this corresponds to Z < Z 2 {5), 
but for the three lowest values of 5, Z 2 {5) < Z < Zi(5) 
and the slope at A = is therefore positive. £l s — f2„ 
nevertheless becomes negative at a finite value of A, and 
in all cases the optimal value of A (for which fl s — f2 n 
has its minimum) is very close to the bulk value Aq. As 
shown by von Delft et al.,i3 only when 5 is very close to 
5 C does the optimal value diminish significantly. It is fur- 
thermore worth noticing that, as can easily be concluded 
from Eq. (^), the optimal value of A is independent of 
Z, for Z < Z x . 

Each time the number of Pauli blocked pairs decreases 
with increasing A, the free energy curve has an upward- 
pointing cusp as exemplified in Fig. ^|by the three small- 
est values of 5. The cusps are spaced more closely at 
smaller values of 5, but the overall envelopes of the three 
curves are rather similar. (The fourth curve has no cusps 
because neither normal nor superconducting states are 
spin-polairzed.) In a model with constant energy level 
spacings, as employed here, the minimum of Q s — f2„ al- 
ways occurs with a minimum of blocked pairs, i.e. in the 
superconducting state the grain will have spin zero for 
even N and spin 1/2 for odd N. 

C. Full T — Z phase diagrams 

At finite temperature both of the functions Z\ (T, 5) 
and Z 2 (T, 5) are found by numerically analyzing Q s — 
ft n as a function of A. In Fig. [| we show the finite 
temperature version of the plots in Fig. |^. The cusps 
are thermally broadened and the condensation energies 
have diminished for all values of Z. The full T-Z phase 
diagram depends on both 6 and electron-number parity. 
Examples at representative intermediate values of 5 are 
presented in FigSr-0 and ^. The general picture is similar 
to the bulk case, even as 5/ A approaches one. At low 



temperature Zi(T,5) > Z 2 (T,5) and the transition to 
the normal state is first order. At the transition A then 
drops abruptly from its Z = value, A(T, 5, Z = 0), to 
zero. At higher temperatures, Z\(T,5) = Z 2 {T,5) and 
the transition is continuous. In that case, the free energy 
curves do not cross zero at finite A if their initial slope 
is positive. 

Whereas Z\(T — > 0, 5) only differs from its bulk value 
in the ultra-small limit, Z 2 (T — > 0, 5) has a significant 5- 
dependence even for intermediate grain sizes. Although 
Z 2 doesn't correspond to a phase transition at low tem- 
peratures, it does have physical meaning as a supercool- 
ing curve, and has heen successfully addressed experi- 
mentally in thin films 2 -!! by measuring enhanced fluctua- 
tions in the neighborhood of Z 2 . 

III. DISCUSSION 

Level spacings near the chemical potential in a real 
ultra-small grain will fluctuated around the mean-value 
used in our idealized model. For an evee, number of elec- 
trons fj,(T = 0, Z = 0) will fall betweenc 2 ! two energy lev- 
els e a and e a _i. rFjpr an odd number (j,(T = 0, Z = Q) 
will fall close tcEa a level e a . Since our results for 
Z 2 (T = 0,5 — * 5 C ) depend primarily on Pauli blocking 
of the first pair of levels, the main consequence of using 
a realistic spectrum are captured by identifying 5 with 
( a - e -i f or cven N and (e a +i - e a -i)/2|Jpr odd N. 
As a consequence, like minimum grain sizes,tj maximum 
Zeeman energies will have a broad distribution. Conse- 
quently, it is difficult to compare directly with specific 
data. Nevertheless our work appears to shed some light 
on the interpretation of both recent and older experi- 
ments. | — | 

In a pioneering early experiment, Giaever and ZellerE2l 
found a violation of the Chandrasekhar-Clogston limit 
which, to our knowledge, is still not fully explained. 
These authors studied tunneling through an ensemble 
of Sn grains with a narrow size distribution. Interpret- 
ing their measurements using a Coulomb blockade pic- 
ture, they concluded that most grains retained a super- 
conducting gap up to magnetic fields that exceeded the 
Chandrasekhar-Clogston limit. The .more recent experi- 
ments by Ralph, Black, and TinkhamB also appear to find 
that the Chandrasekhar-Clogston limit can be exceeded. 
In the later tunneling experiments, a quasiparticle gap, 
is observed to decrease linearly with the Zeeman cou- 
pling strength Z . The linear dependence arises from the 
Zeeman splitting of quasiparticle energies and is consis-. 
tent with the mean-field theory employed in this paper .til 
The linear decrease continues to Zeeman fields that ex- 
ceed the Chandrasekhar-Clogston limit without the dis- 
continuous drop which would be expected if A dropped 
abruptly to zero. However, this observation is made at 
a value of Aq/5 which is smaller than those for which 
the Chandrasekhar-Clogston limit is exceeded in a model 
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with equally spaced energy levels. The experiment could 
be explained by assumming that this particular sample 
happens to have a relatively large energy spacing at the 
Fermi energy. 

The results presented in this paper are based on mean- 
field theory, and a few cautionary remarks concerning 
its validity are in order. For T = and Z < Z%, the 
mean-field condensation energy for 8 — > 8 C becomes mi- 
croscopic: Q s — fi„ — > — S[ln(S c /8} 2 /((3, a). Thermal 
fluctuations will therefore be important unless the tem- 
perature is well below— the bulk critical temperature and 
quantum fluctuations^ will be increasingly important as 
8 C is approached. Clearly the ultrasmall grain in this 
regime will not exhibit anything approaching a true phase 
transition between normal and superconducting states. 
The phase boundaries found in mean-field-theory should 
be regarded as estimates for the localtions of crossovers 
which become more gradual as 8 increases. 
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A. D. Zaikin are gratefully acknowledged. 



1 A critical examination of the assumption that pairing oc- 
curs only between identical orbitals will be presented else- 
where. M. C. B0nsager and A. H. MacDonald, in prepara- 
tion, to be submitted to Physical Review B. 

17 A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962); B. S. Chan- 
drasekhar, Appl. Phys. Lett. 1, 7 (1962). 

18 P.-G. de Gennes and M. Tinkham, Physica 1, 107 (1964). 

19 A. L. Fetter and J. D. Walecka, Quantum Theory of Many- 
Particle Systems, McGraw-Hill, New York, 1971. 

20 G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963); K. Maki, 
and T. Tsuneto, Prog. Theor. Phys. 31, 945 (1964). 

21 R. Meservey and P. M. Tedrow, Phys. Rep. 238, 173 
(1994). 

22 For a constant energy level spectrum /i is independent of 
T and Z and falls exactly on a level (odd N) or exactly 
between two levels (even N). This is not the case for a 
random spectrum. Rather, fi will depend on both T and 
Z and it will not fall exactly on or between levels. M.C. 
B0nsager and A.H. MacDonald, unpublished. 

23 I. Giaever and H. R. Zeller, Phys. Rev. Lett. 20, 1504 

(1968) ; H. R. Zeller and I. Giaever, Phys. Rev. 181, 789 

(1969) . 

24 K. A. Matveev and A. I. Larkin, Phys. Rev. Lett. 78, 3749 
(1997); A. Mastellone, G. Falci, and R. Fazio, Phys. Rev. 
Lett. 80, 4542 (1998); S. D. Berger and B. I. Halperin, 
Phys. Rev. B 58, 5213 (1998). 



2 



1 M. T. Tuominen, J. M. Hergenrother, T. S. Tighe, and M. 
Tinkham, Phys. Rev. Lett. 69, 1997 (1992). 

2 P. Lafarge, P. Joyez, D. Esteve, C. Urbina, and M. H. De- 
voret, Phys. Rev. Lett. 70, 994 (1993). 

3 T. M. Eiles, J. M. Martinis, and M. H. Devoret, Phys. Rev. 
Lett. 70, 1862 (1993). 

4 D. V.Averin and Yu. V. Nazarov, Phys. Rev. Lett. 69. 1993 
(1992) 

J B. Janko, A. Smith, and V. Ambegaokar, Phys. Rev. B, 
50, 1152 (1994). 

6 D. C. Ralph, C. T. Black, and M. Tinkham, Phys. Rev.. 
Lett. 74, 3241 (1995). 

7 C. T. Black, D. C. Ralph, and M. Tinkham, Phys. Rev. 
Lett. 76, 688 (1996). 

8 D. C. Ralph, C. T. Black, and M. Tinkham, Phys. Rev. 
Lett. 78, 4087 (1997). 

9 P. W. Anderson, J. Phys. Chem. Solids 11, 26 (1959). 

10 O. Agam, N. S. Wingreen, B. L. Altshuler, D. C. Ralph, 
and M. Tinkham, Phys. Rev. Lett. 78, 1956 (1997). 

11 O. Agam and I. Aleiner, Phys. Rev. B 56, R5759 (1997). 

12 J. von Delft, A. D. Zaikin, D. S. Golubev, and W. Tichy, 
Phys. Rev. Lett. 77, 3189 (1996). 

13 R. A. Smith and V. Ambegaokar, Phys. Rev. Lett. 77, 4962 
(1996). 

14 F. Braun, J. von Delft, D. C. Ralph, and M. Tinkham, 
Phys. Rev. Lett. 79 , 921 (1997). F. Braun and J. von Delft, 



cond-mat/9801170 



R. Balian, H. Flocard, and M. Veneroni, 



and cond-mat/9802006 



nuc-th/9706041 



— Z,(T=0) 
- Z 2 (T=0) 
--- approx. Z,(T=0) 



< 1 
N 



0-5 i'fl/ 



' / 



5/1 



FIG. 1. The zero temperature limit of Z\ (solid line) and 
Zi (long-dashed line) as a function of <5 for an even number of 
electrons. The dashed line is the approximate expression (^) 
for Z\ derived in the text. Z-2 can be evaluated analytically 
as explained in the text. The discontinuities in this curve 
occur at 8/A = 0.4821, 0.2475, 0.1659, 0.1247, ... and the 
corresponding minima and maxima are given by Eq. (^) . The 
function Z\(T — 0,S) is continous but it has discontinuities 
in its first derivative at a series of 8. This fact is further 
illustrated in Fig. ^. 
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FIG. 2. As in Fig. |l| but for 

ber of electrons. The discontinuities 
5/A = 0.3276, 0.1987, 0.1424, 0.1109, .... As it is the case 
for even N (see Fig. [fj), the function Zi(T — 0,5) has dis- 
continuities in its first derivative at a series of 5 given by 
Zi(5) = (n — a)S, n = 2,3,4, i.e. where extrapolations 
of Z2 (indicated by dotted lines) cross Z\. These cusps are 
a consequence of discontinuities of the first derivative of the 
normal state energy with respect to Z. In most cases, the 
cusps are too small see in this figure, however. 
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The same as in Fig. ^ but at finite temperature 
0.2, where T c o is the bulk critical temperature at 
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FIG. 3. The dependence of the free energy difference per 
particle (fi s - fi„)/iV on A at (T = 0, Z = 0.65A ) for four 
different values of level spacing 5 for an even number of elec- 
trons. The circles indicates cusps in the curves. The case 
<5/Ao = 1.4 is special in the sense that it corresponds to Zee- 
man field below Z 2 (5) (see Fig. §). For Z > Z a (1.4A ) this 
curve would also have a positive slope at A = 0. 
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FIG. 5. The functions Zi(T) and Z 2 (T) for both even 
and odd numbers of electrons. This figure is for 5/Aq = 0.4. 
Temperatures are expressed in terms of the bulk critical tem- 
perature Tco. For S/Aq = 0.4 Z\(T = 0)/Aq is already close 
to its bulk value, l/y/2. For both even and odd number of 
electrons the transition becomes first order near T/T c o = 0.56, 
close to the corresponding bulk value. Z2 shows the largest 
deviation from the bulk limit. For odd number of particles 
Z 2 {T = 0,5) < Z 2 (T = 0,S -> 0) while for even number of 
particles Z 2 (T = 0, 5) > Z 2 (T = 0, 5 -> 0). 
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FIG. 6. The functions Zi(T) and Z 2 (T) for both even 
and odd numbers of electrons. This figure is for 5/Ao = 0.79, 
a value sufficiently large to yield phase diagrams which dif- 
fer substantially from their bulk counterparts. The transi- 
tion becomes first order at T/T c o = 0.76 for even N and at 
T/Tco = 0.44 for odd N. This illustration also reflects the de- 
pendence of the critical temperature 8 and parity discussed by 
by von Delft et al.13. For odd N T c decreases monotonically 
with 8 whereas for even N it increases up to about S/Ao ~ 2 
before it decreases. 
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